


•LOAN copy: return to 
AFWL TECHNICAL LIBRARY 
KIRTLAND AFB,; N. M. 


FIXED-RANGE OPTIMUM TRAJECTORIES 
FOR SHORT-HAUL AIRCRAFT 

Heinz Erzberger, John D. McLean, 
and John P. Barman 

Ames Research Center 
Moffett Field, Calif. 94033 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • DECEMBER 1975 



TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 



0133*133 


1. Report No. 2. Government Accession No. 

NASA TN D-8115 

3. Recipient's Catalog No. 

4. Title and Subtitle 

FIXED-RANGE OPTIMUM TRAJECTORIES FOR SHORT-HAUL AIRCRAFT 

5. Report Date 
December 1975 

6. Performing Organization Code 

7. Author(s) 

Heinz Erzberger, John D. McLean, and John F. Barman 

8. Performing Organization Report No. 
A-5874 

10. Work Unit No. 
505-07-11 

9. Performing Organization Name and Address 
NASA Ames Research Center 
Moffett Field, Calif., 94035 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Note 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 

John F, Barman is an NRC Research Associate at Ames Research Center 


16, Abstract 

An algorithm, based on the energy-state method, is derived for calculating optimum 
trajectories with a range constraint. The basis of the algorithm is the assumption that optimum 
trajectories consist of, at most, three segments: an increasing energy segment (climb); a constant 
energy segment (cruise); and a decreasing energy segment (descent). This assumption.. allows energy 
to be used as the independent variable in the increasing.;and decreasing energy segments, thereby 
eliminating the integration of a separate adjoint differential equation and simplifying the 
calculus of variations problem to one requiring only pointwise extremization of algebraic functions. 
The algorithm is used to compute minimum fuel, minimum time, and minimum direct-operating-cost 
trajectories, with range as a parameter, for an in-service CTOL aircraft and for an advanced STOL 
aircraft. For the CTOL aircraft and the minimum-fuel performance function, the optimum controls, 
consisting of air-speed and engine power setting, are continuous functions of the energy in both 
climb and descent as well as near the maximum or cruise energy. This is also true for the STOL 
aircraft except in the descent where at one energy level a nearly constant energy dive segment 
occurs, yielding a discontinuity in the airspeed at that energy. The reason for this segment 
appears to be the relatively high fuel flow at idle power of the engines used by this STOL aircraft. 
Use of a simplified trajectory which eliminates the dive increases the fuel consumption of the 
total descent trajectory by about 10 percent and the time to fly the descent by about 19 percent 
compared to the optimum. 


17. Key Words (Suggested by Author(s) ) 


18. Distribution Statement 


Minimum fuel trajectories for aircraft 


Shott-haul aircraft 
Energy state method 


STAR Category — 

Unlimited 

04 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 


31 

$4.25 


For sale by the National Technical Information Service, Springfield, Virginia 21161 




SYMBOLS 


C 



Cip 

D(L,h,V) 

E 

E 


h 

I 

J 

k(m) 

k 

L 

M 

P 


direct operating cost (DOC), dollars/trip 

constant cost component of DOC, dollars 

fuel-independent component of DOC, dollars 

drag coefficient 

zero-lift drag coefficient 

unit cost of fuel, cents /kg 

lift coefficient 

cost per unit flight time, dollars/hr 
drag, N 

total energy, m 
energy rate, m/sec 
fuel consumed, kg 

normalized fuel flow function, kg/hr 
normalized thrust force function, N 
acceleration of gravity, m/sec^ 
altitude, km 
integrand 

integral performance index 
induced drag factor 
lift curve intercept 
slope of the lift curve 
lift, N 
Mach number 

integrand performance index 


iii 



p 


R 

R 


dn 




ndn 


nup 


R 

R 

r 

S 

T 

t 

V 

W 
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fuel tradeoff coefficient 


T atmospheric temperature, °K 

standard sea level atmospheric temperature, °K 
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f final 

i initial 

in inlet 
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opt optimum 
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FIXED-RAUGE OPTIMUM TRAJECTORIES FOR SHORT-HAUL AIRCRAFT 


Heinz ErzUerger, John D. McLean, and John F. Barman^ 
Ames Research Center 


SUMMARY 


An algorithm, based on the energy- state method, is derived for calculat- 
ing optimum trajectories with a range constraint. The basis of the algorithm 
is the assumption that optimum trajectories consist of, at most, three seg- 
ments: an increasing energy segment (climb); a constant energy segment 

(cruise); and a decreasing energy segment (descent). This assumption allows 
energy to be used as the independent variable in the increasing and decreasing 
energy segments, thereby eliminating the integration of a separate adjoint 
differential equation and simplifying the calculus of variations problem to 
one requiring only pointwise extremization of algebraic functions. The algo- 
rithm is used to compute minimum fuel, minimum time, and minimum direct- 
operating-cost trajectories, with range as a parameter, for an in-service CTOL 
aircraft and for an advanced STOL aircraft. For the CTOL aircraft and the 
minimum- fuel performance function, the optimum controls, consisting of air- 
speed and engine power setting, are continuous functions of the energy in both 
climb and descent as well as near the maximum or cruise energy. This is also 
true for the STOL aircraft except in the descent where at one energy level a 
nearly constant energy dive segment occurs, yielding a discontinuity in the 
airspeed at that energy. The reason for this segment appears to be the rela- 
tively high fuel flow at idle power of the engines used by this STOL aircraft."^ 
Use of a simplified trajectoiy which eliminates the dive increases the fuel 
consumption of the total descent trajectory by about 10 percent and the time 
to fly the descent by about 19 percent compared to the optimum. 


INTRODUCTION 


Sharply escalating fuel prices and the threat of future fuel shortages 
have generated strong interest in finding methods of reducing aviation fuel 
consumption that do not seriously affect the level of airline service. One 
such method, flight -path optimization, has the potential for saving signifi- 
cant quantities of fuel. Although there is a long history of flight-path 
optimization studies for all types of aircraft, modern approaches based on the 
variational calculus have been applied more frequently to supersonic military 
rather than subsonic civil aircraft missions. 

Subsonic aircraft missions can be roughly divided into long haul and 
short haul. In long-haul missions, which are characterized by long periods of 
cruise flight, the central problem in flight planning is optimizing the ground 
track and the altitude profile during cruise so as to use wind, temperature. 
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and other atmospheric conditions to the greatest advantage. For flying the 
climbout and descent segments , long-haul operators use procedures supplied by 
the aircraft manufacturers. Such procedures are generally not optimum, in 
that they do not minimize a performance index such as time or fuel used and, 
because they affect only a small portion of the total flight path, they have a 
limited in 5 )act on long-haul flight performance. 

These priorities are essentially reversed in short-haul missions (800 km 
or less). For these missions, the cruise segment is relatively short and 
therefore the climbout and descent segments play the dominant role in flight- 
path optimization. Moreover, the short range of these missions can yield a 
high degree of interdependence of the climbout, descent, and cruise optimiza- 
tion problems, suggesting that range should enter explicitly as a boundary 
value . 


This report describes a concept ion ally and computationally simple algo- 
rithm for calculating optimum flight paths with a range constraint. The algo- 
rithm is applied to flight profile optimization of two t 3 rpes of short-haiil 
aircraft: a currently in-service CTOL jet and a future-design jet STOL air- 

craft (an augmentor wing with a supercritical airfoil). Minimum fuel and 
time-flight paths are computed for each aircraft. Furthermore, assuming that 
the relative costs of time and fuel are known, a procedure is given for 
selecting minimum direct-operating-cost (DOC) flight paths. 


DEFINITION OF PERFORMANCE INDEX 


For a fixed aircraft configuration and a specified origin and destina^- 
tion, the primary factor that influences aircraft performance is operating 
procedure. For an airline, the performance goal is generally minimization of 
the direct operating cost (DOC) which consists of fuel cost, crew cost, main- 
tenance cost, depreciation, and insurance. It is therefore of interest to 
detemiine a relationship specifying the dependence of DOC on flying time and 
fuel consumed. Relating DOC to fuel consumed is straightforward, but relating 
the other components of DOC to flying time is more difficult. A reasonable 
approach for the purpose of developing optimum flight procedures is to param- 
eterize the cost of the DOC components, other than fuel, by the relation 

Cl = Co + C^T (1) 

where C^ is a cost component which is independent of time over the time 
scale of the mission, C«p is the cost per unit flight time, and T is the 
total flight time (ref. l). Thus, the factor C«p should reflect the monetary 
value of flight time, exclusive of fuel cost. The total cost of the flight 
can therefore be written as 


C = CpP + C^T + Co 

or in integral form as 

C = J (CpWp + C^)dt + c^ 


( 2 ) 

(3) 


2 



where Cp is the unit cost of fuel and Wp is the fuel flow rate. Since min- 
imization of equation (3) depends only on the ratio Cp/Cij and not on an 

equivalent performance index, which is more convenient for performance 
optimization, is 


where 


J = [ Wp + (l - a)]dt ; 0 < a < 1 

Cp/ Cip 

= 1 + {c^/c^y 


(U) 

(5) 


The two extreme values of o give the two important special cases of minimum 
time (a = O) and minimum fuel (a = l) performance indices. 


DERIVATION OF OPTIMIZATION ALGORITHM 


The approach to flight-path optimization talcen here follows the trend in 
the recent literature (refs. 2 and 3) of using the energy state formulation as 
the basis for computing the optimum flight paths. The rate of change of 
energy along the flight path is given by 

H = [Tp(n,h,v) - D(L,h,v)] I (6) 


where 


h = E 



2 


and Tp is the thrust, D the drag, L the lift, W the weight of the air- 
craft, h the altitude, V the airspeed, and g the acceleration of gravity. 
The controls which determine the flight path are taken as the power setting, 
n, and airspeed V. In accordance with the energy-state method, it will be 
assumed that lift is equal to weight in computing the drag. The flight path 
is also constrained to cover a specified range, R, which makes it necessary to 
introduce an additional state equation: 


dR 

dt 


V 


(7) 


In equation (?) it is assumed that, except at isolated time instants, the 
flight-path angle, y, is small - allowing the small-angle approximation, 
cos Y = 1, to be made - and that the wind speed relative to the airspeed is 
negligible. The latter assumption is not necessary for the development that 
follows, but the effect of winds on the optimum flight path is not considered 
in this report. Finally, the energy and range are specified at the beginning 
and end of the flight path: 
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E = Ei 

and 

II 

at 

t = 0 

E = E^ 

and 

II 

at 

t = T 


( 8 ) 


Minimization of equation (U) subject to the constraints of equations (6) 
and ( 7 ) and the boundary conditions of equation (8) is a frequently studied 
problem in optimum control. However, solutions published in the literature 
apply mostly to military supersonic aircraft missions. The most recent 
results in the field have generalized the energy-state method to include 
flight with turns, again for military aircraft applications. References 2, 3, 
and 4 provide a brief list of recent papers on this subject. 

Application of the maximum principle results in a fourth-order, two-point 
boundary value problem consisting of the two state equations (6) and (7), the 
associated adjoint equations, and the boundary conditions (8). Numerical solu- 
tion of two-point boundary value problems of this type has proven to be diffi- 
cult. Zagalski (ref. 4) described a procedure for simplifying the two-point 
boundary value problem in the energy-state formulation. In this report, an 
assumption about the structure of the optimum flight paths is similarly used 
to eliminate the need for the adjoint equation. The assumption is as follows: 

Optimum flight paths for the problem defined here consist of three seg- 
ments at most; in the first, energy increases monotonically , in the second it 
is constant along with the velocity, and in the third it decreases monotoni- 
cally with time. The resiilting flight paths may not be optimum for the prob- 
lem as originally stated but should be very nearly so for most problems of 
interest. 


To show how the assumption simplifies the confutation of optimum flight 
paths, equation (4) is written as the sum of costs of the three segments: 


up 

J = \ * (P)4^„(T^-T^p) * 


1 




(9) 


dn 


where P = aWn> + (l-a), T is the time at the end of the increasing energy 
^ up 

segment, and T, 3 ji is the time at the start of the decreasing energy segment. 
The rate of change of energy, £, is given by equation (6). Since airspeed is 
constant in cruise, the middle term corresponding to the cruise cost can be 
expressed as 


^^^E=0^^dn “ 


T ) = 
up 


XR 


( 10 ) 


In equation (lO), R^, is the distance spent in cruise and A, defined as the 
cruise efficiency, is evaluated at the cruise energy using the relationship 

It 





where is the cruise airspeed. The cruise distance 


is given by 


where 



( 12 ) 


(13) 


and Rq must satisfy the inequality R^ > 0. Since by the assumption of the 
preceding paragraph, energy changes monotonically in the first and last terms" 
of equation (9)» energy can replace time as the independent variable in equa- 
tions (9) and (13), After applying the transformation dt = dE/tl to these 
equations and substituting the appropriate integration limits, they become 


J 


R 

un 




(lU) 


(15) 


Here, ^refers to the maximum energy (cruise takes place at Ejj^ax)* Note 

from equation (lU) that since R^ > 0, the cruise efficiency X should be 
chosen as small as possible in the process of minimizing J described below. 
Equations (l2), (l4), and (l5) can be combined into a single equation: 


i E 

max / V 

dE + AR (16) 

If we assume for the moment that known, then equation (l6) can 

be minimized by performing three independent algebraic minimizations. The 
first optimizes cruise conditions and consists of minimizing X at E^^^x 
not^d above. The remaining two problems are the minimization of the two inte- 
gral cost terms in equation (l6). They are minimized by choosing power set- 
tings and airspeeds as functions of E so that each of the two integrands is 
minimized at all values of the independent variable E throughout its inte- 
gration interval. The minimization of each integrand must obey different con- 
straints on the controls, namely, II and V must be such that > 0 for the 
increasing energy segment and < 0 for the decreasing one. In addition, 
there are constraints on the power setting and airspeed which depend on the 


J = 


! E 

max 


(P-XV) 


E. 

1 



energy. These minimizations , defined below, ensure that J is minimum for 

each choice of E : 
max 


X . (e ) 

opt max 


I 

up 


(E,E ) 
max 


min 

n,v 


min 

n,v 




E fixed 


I, (E,E ) 
dn ’ max 


min 

n,v 


A<o 

E fixed 


(17) 


( 18 ) 


(19) 


The last step in minimizing J involves finding the optimum E^^^x 
the specified range. A necessary condition for J to be minimum at some 

Is foT the dertvativo of J with respect to Ejj^ax vanish at Ejj^ax* 
This is done as follows. We compute the derivative of J for an increasing 
sequence of E^^^y , starting with E^^ ^ = max{Ej_,E^}. As Ejj^ax increased 
for fixed Ej Rq decreases and two cases can occur: 

(1) 6J = 0 for some Ejj^ax s'^ch that 0 < < R or 6J 0 as R^ 0. 

(2) 6J < 0 as ^ 0. 

We shall discuss case (l) in more detail than case (2) since case (l) applies 
to both aircraft models studied in the following sections of this report. 


Case (l) 


For this case, 
puted from equation 
( 17 ) through ( 19 ) : 


the derivative of J with respect to com- 

(lU) using Leibnitz’s rule and equations tl2), (l5), and 


dJ 

dE 

max 


= l“ (E ,E ) + I" (E ,E ) + R 
up max max dn max max c 



dE 


E=E 

max 


( 20 ) 


The quantities I and are limiting values of I and I. as E E 

^ up dn ° up dn max 

from below. In general, the limits must be evaluated for this case since both 

the numerators and denominators of the bracketed quantities in equations (l8) 

and ( 19 ) can approach zero simultaneously as E ^ E 

max 

For the two types of aircraft studied in this report, the fxmctions ^up» 
and have special characteristics that have been observed in 
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computations and have also been verified analytically using polynomial repre- 
sentations for the aerodynamic and engine models given in the next section. 
The derivation of some of these characteristics is given in the appendix. 

First, the function X . (E ) has the general shape shown in figure 1, 

’ opt max 

Second 5 the functions and obey the relationship 


1“ (E ,E ) + (E ,E ) = 0 
up max’ max dn max’ max 


(21a) 


^opt 



Thus the values of the limit func- 
tions are equal in magnitude and 
opposite in sign. Moreover, at 
Ejjj^ax optimum controls satisfy 

the relationships: 


V 


opt 

opt 


= Arg. min — 

n,v 


E=0 


p - 


Figure 1,- Typical relationship 
between cruise efficiency and 
cruise energy; subsonic aircraft. 


= Arg. min lim 

= Arg, min lim 


X .V 
out 


E 


P - 


X . V 
opt 


E 


E>0 


E<0 


This relationship is illustrated by figure 2 in which is plotted the 
optimum controls obtained from equations (l8) and (l9) as a function 


(211d) 

locus of 
of E. 


P-XV 


7T 


MAX- 


POWER 

SETTING 

OR 

THRUST 




'MIN 





V E )-XV=0 

'''max' 

PERMISSIBLE 

CONTROL 


/////7777777Z^///////////////A region 


'MIN 


MAX 


Vc 

AIRSPEED 

LOCI OF OPTIMUM CONTROLS CONVERGE AT 


E^ax AIRCRAFT STUDIED 


Figure 2.- Typical locus of optimum controls as a function of E, 
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It is seen that the two hranches of the locus converge to the controls 
obtained from equation (l7) for cruise at Emax* Furthermore, it can he shown 
by analysis of the polynomial representation that the bracketed quantities in 
equations (l8) and (ip) have a local minimum at with respect to admis- 

sible control variations. However, the continuity of the optimum controls at 
Emflv depends on the aerodynamic characteristics and especially on the rela- 
tionship between thrrist and fuel flow as derived in the appendix. In particu- 
lar, it is not generally true for the linear dependence often assumed in 
earlier work (refs. 3 and U). 


Substituting equation (2l) into (20 ) yields the sinplified form 

dX 


dJ 


dE 


= R 


max 


opt 

dE 


= 0 


E=E 


max 


( 22 ) 


It follows immediately from equation (22) that a nonzero crtiise distance can 
be optimum only if Emay yields the minimtim of Xq ^(e), which is designated 
as X* in figure 1. This resuilt is consistent with the results of other 
workers, for example, Schultz and Zagalsky (ref. 3). If the range is less 
than some minimum value, the maximum energy of the flight path will fall below 
the optimum cruise energy. In that case, equation (22) requires zero cruise 
distance (R^ = O). This fact simplifies computing the relationship between R 
and Eniay . For each choice of E^^y , one first conputes ^opt^^ux) sind then 
Rup and using relations (15). In performing the range integrations, the 

optimum controls obtained from equations (l8) and (l9) must be used. The 
range-energy relationship is then 


r(e ) 

= R 

+ R. . 

E 

< E ^ 

' max 

up 

dn * 

max 

optc 

r(e ) 

= R 

+ R^ + R , 

E 

= E ^ 

max 

up 

dn c * 

max 

optc 


I 


(23) 


By integrating equations (l8) and (ip) for each choice of E^n^y , the relation- 
ship between optimum cost and E^ay can also be obtained: 


J *(E ) = 

opt max 


) .E 

max 

1 


i E 

max 

E„ 




(E 


max' optc' max) 


(2l^) 


Case (2) 


In this case the quantity ^op-t loses its previous interpretation as a 
measTxre of cruise efficiency and instead becomes an adjoint variable X, whose 
value mTist be chosen iteratively to achieve a specified range R. The condi- 
tion which determines the maximiun energy for a trial value of X can be 
obtained by application of the maximum principle to equations (lU) and (15) 
with Rq set to zero (ref. 5)» The condition is found to be 
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min 

n.v 


0 


/e>0 


E=E 

max 



E=E 


max 


For this case the optimum controls will generally not be continuous at the 
maximum energy point of the trajectory. 

An essential requirement for implementing the algorithm on a computer is 
a routine for minimizing a function of two variables. Results presented in 
the next sections were obtained using an algorithm known, in the literature af 
nonlinear programming, as the method of local variations (ref, 6), This opti- 
mization algorithm does not require derivative calculations. The minimization 
over power setting can be eliminated a priori only for the minimum time per- 
formance index, a = 0, since in that case the power setting is a bang-bang 
function of energy. Special provisions must also be made for computing I^p 
and as E approaches Ejj^ax where the optimum controls yields values of £ 

near zero. A satisfactory approach is to decrease the integration step size 
as E Ejj^ax stop when |£l < Ernin* ^ value of 0,33 m/sec for E^iin 

was used in the results described in the following section. 


AERODYNAMIC AND PROPULSION MODELS 


The algorithm was used to study optimum trajectories for two types of 
aircraft: the CTOL aircraft, a l80-passenger tri-jet currently in short-haul 

service; and the STOL aircraft, a future-design, four-engine, 150-passenger, 
TOO-m field length jet aircraft using augmentor flaps to achieve STOL perfor- 
mance, and a swept, supercritical airfoil to achieve a Mach 0,8 cruise speed. 
The design of the STOL was carried out in a recent NASA study (ref. 7). 

Aerodynamic forces for both types of aircraft are described by the 
following equations: 


L = Cj^Sq 

(lift force) 

(25) 

D = Cj^Sq 

(drag force) 

( 26 ) 


(lift coefficient) 

(27) 


\ 

where Cj) is the drag coefficient, S the wing reference area, and q the 
dynamic pressure; q = (l/2)p(h)V^, The air density p , as a function of alti- 
tude h, was obtained from the 1962 Standard Atmosphere. The drag coeffi- 
cients for the two aircraft are parameterized as follows: 
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.2 

"L 


0.1)2 

(CTOL aircraft) 

(28) 

a(m,Cl) 

(STOL aircraft) 

(29) 


where M is Mach mimher. 


The constants in equation (27) and the Mach number corrections in equa- 
tions (28) and (29) together with the masses and wing reference areas for the 
two aircraft are given in tables 1 and 2. 


TABLE 1.- AERODYNAMIC DATA AND WEIGHT FOR CTOL AIRCRAFT 


M 

0 

' O.TO' 

‘ 0 . 7 '^ 

0.80 

0.82 

■ 0.8¥“' 

0.85 

0.86 

"0.87 

0.88 

0.89 

Cd^(M) 

K(M) 

0.0173 

0 . 0861 t 

0.0173 

O.O86U 

0.0174 

0.0932 

O.OITT 

0.103 

0.0179 

0.113 

0.0182 

0.128 

0 .' 0 l 86 

0 .l 4 i 

0.0189 

0.169 

0.0195 

0.184 , 

0.0203 

0.211 

0.0218 

0.241 


Note: ICq = 0, ki = 0.1, S = IU5 and VJ = 68,200 kg. 


TABLE 2.- AERODYNAMIC DATA AND MASS FOR STOL AIRCRAFT 




0.62 

0.64 

0.66 

0.68 

0.70 

0.72 

0.74 

0.76 

0.77 

Mach 

numher and lift 

coefficient correction to drag coefficient, A(M,Cl) 

0 

0 

0.0001 

0.0002 

0.0005 

0.00075 

0.0010 

0.0013 

0.0017 

0.0026 

0.0034 

.20 


.0001 

.0002 

.0005 

.00075 

.0010 

0.0013 

.0017 

.0026 

.0034 

.25 



.0001 

.0002 

.0005 

.0009 

0.0012 

.0016 

.0021 

.0026 

.30 




.0001 

.0003 

.0007 

0.0010 

.0014 

.0020 

.0025 

.35 




.0001 

.0003 

.0007 

0.0011 

.0014 

.0020 

.0024 

.ho 




.0001 

.0003 

.0007 

0.0011 

.0015 

.0021 

.0025 


Note: = -0.1, k^ = 0,075, S = 255 , and W = 90,909 kg. 


The aircraft propulsion systems are modeled using corrected engine 
parameters (ref. 8) as follows: 

T 

•^ = f^(n//e,M) (30) 

— = f (n//e,M) (31) 

6/0 ^ 

where Tp and Wp are the thrust and fuel flow rate, respectively. The quan- 
tities 6 and 0 are pressure and temperature factors given by 
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6 


( 32 ) 


0 


p(h) 



Y-1 

2 


y/(y-i ) 


(It m2) 


(33) 


where x(h)/TQ and p(h)/pQ are atmospheric temperature and atmospheric pres- 
sure ratios obtained from the 1962 Standard Atmosphere and y = 1,U is the air 
specific heat ratio. The q^uantity II is the power setting, which for the CTOL 
aircraft is actual rotor revolution per minute and for the STOL aircraft is 
defined as 


n = /e (3U) 

In 

where Tj^ is the actual turbine inlet temperature. The corrected thrust and 
fuel flow curves per engine for the CTOL and STOL aircraft are given in fig- 
ures 3 and U, respectively. (in figure and in subsequent discussions of 
STOL trajectories, II is referred to as T^ , the corrected turbine inlet 
temperature. ) 


Fixed-Range Optimum Flight Paths: CTOL Aircraft 

The first step in the implementation of the algorithm is the computation 
of cruise efficiency as a function of maximum energy (eq. (iT)) for different 
values of a. The results of this computation for a fixed mass of 68,200 kg 
are plotted in figure 5 in cruise-fuel-efficiency (kg/km) and cruise-airspeed 
(m/sec) coordinates. The solid heavy line gives the envelope of optimum 
cruise conditions as a ranges from zero to 1, Points A and B define the 
fuel efficiencies and airspeeds for minimum fuel and minimum time, respec- 
tively. Minimum- fuel cruise occurs at an altitude of 10 km, which is about 
1 km below the ceiling, whereas minimum-time cruise (maximum airspeed) occurs 
at an altitude of 5.5 km. Figure 5 also shows the loci of fuel efficiency and 
airspeed for o = 0.5 as the maximum energy is allowed to vary over its pos- 
sible range. It was shown earlier that no cruise can take place in the opti- 
mum trajectories except at the absolute optimum cruise represented by points 
on the solid line. Point C on the a = 0.5 cruise locus achieves this optimum 
value denoted by X*. Any other point on the a = 0.5 cruise locus corresponds 
to a ^opt value greater than X* and therefore corresponds to a range less 
than the smallest range containing a nonzero cruise distance. 

In computing the cruise efficiency, the actual rotor revolutions per 
minute were restricted to maximum cruise power setting of the engine which is 
modeled as : 


n 

max. 

cruise 


11 


7.48 + 0.l8(h/3 - 1.5) , h < 7.62 km 
7.66 otherwise 


( 35 ) 



3456789 10 

Tr/ye, krpm 


(a) Corrected thrust as a function of Mach number 
and corrected pover setting. 



23456789 


ir/ye^ krpm 

(b) Corrected fuel flow rate as a function of Mach 
number and corrected power setting. 

Figure 3.- CTOL aircraft: thrust and fuel flow rate corrections 



CORRECTED THRUST Fn/8, NxIO" 


I 





CORRECTED TURBINE INLET TEMPERATURE 



CORRECTED TURBINE INLET TEMPERATURE (T|/0), «K 


(a) Corrected thrust as a function of (h) Corrected fuel flow rate as a 
corrected turbine inlet temperature. function of corrected turbine 

inlet temperature. 


Figure U.- STOL aircraft: 


thrust and fuel flov rate corrections. 



Figure Cruise operating cost efficiency, CTOL aircraft. 
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In computing the optimum climb and descent strategies, dynamic pressure 
q is limited to 35 N/m^, the altitude is constrained to be larger than zero, 
and the actual engine power setting is bounded from below by 3*95 krpm* For 
0,2 < 0 ^ 1 , the power setting is limited to maximum climb which is modeled as 


n T.TO + 0.02Th , h > 7.62 km 

max. = - 

climb 7.90 otherwise 


(36) 


For 0 < 0 ^ 0.2, which corresponds to near minimum-time missions, the power 
setting is limited to maximum cruise during the entire flight profile. This 
limitation will exclude prolonged operation of the engine at maximum allowable 
power setting which, from the point of view of cost of maintenance, might be 
undesirable. 


The climb and descent range, time, and fuel are calculated by integrating 
equations (l5) and the following two equations using the trapezoidal rule: 


T 

up 
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(37) 


(38) 


It was shown earlier that as E ^ax» optimum power setting and 

optimum velocity approach the cruise values corresponding to ^opt^^ax^s that 
is, E approaches zero causing the integrands of equations (l5), (37), and (38) 
to diverge. To circumvent this problem, a decreasing step size for AE is 
employed as E ^ ^ax» criterion used for obtaining an appropriate inte- 

gration step size is given by 


At = < 30 sec (39) 

|E| 

where At is the time increment. The integration process is stopped when 
|fi| < 0.3 m/sec. The sensitivity of the minimum-fuel trajectories with 
respect to At was studied. If At is restricted to be less than 15 sec, 
rather than 30 sec, then the total time and fuel of the mission for a fixed 
range are changed by only 0.3 percent. 

Figure 6 shows optimum climb and descent trajectories in altitude-velocity 
coordinates for minimum fuel ( 0 = 1 ), minimum time ( 0 = 0 ), and an intermediate 
value of 0 = 0.5. For each of the three values of 0 , optimum trajectories 
corresponding to different total ranges are plotted, with the range indicated 
on each trajectory. For each value of 0 , the trajectory with the largest 
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ALTITUDE, km 


cr=| cr=.5 



Figure 6,- Speed-altitude profiles, 
CTOL aircraft. 


value of Ejj^ax optimum 

trajectory corresponding to the 
smallest range with a nonzero 
cruise segment. The only pos- 
sible cruise points are also 
shown in figure 6. Note that 
the velocity altitude profile 
^max is ^ i nuous for 
each optimum trajectory. 

Figure 7 plots altitude 
versus range of several opti- 
mum trajectories 5 each for a 
different a value. All tra- 
jectories cover a range of 800 
km. As a approaches zero, 
the optimum cruise altitude 
decreases while the range 
covered in the climb trajectory 
increases . 



Figure 7.- Altitude-range profiles, 
CTOL aircraft. 


Normalized throttle vari- 
ations versus normalized range 
are plotted in figure 8 where 
n is the actual rotor revolu- 
tions in krpm; = 3.95 

krpm is the idle throttle set- 
ting; is maximum climb 

throttle setting for 0.2 a ^ 

1 and is the maximiom cruise 
throttle setting for 

0<a<0.2;R is the dis- 

nup 

tance from the maximum energy 

point divided by R and R _ 

^ up ndn 

is the equivalent fraction of 
R(3^. The intersection of these 
curves with the throttle set- 
ting axis corresponds to the 
cruise throttle setting. It 
is seen that as a ^ 0 the 
fraction of the climb range 
covered at maximum throttle 
setting increases and the frac- 
tion of descent range covered 
at idle throttle setting 
decreases. This is expected 
as time is weighed more and 
fuel less in the performance 
index. The portion of the 
optimum trajectories 
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Figure 8.- Percent throttle setting versus percent climb or 
descent range covered vith an 800-km fixed range; CTOL 
aircraft. 

corresponding to intermediate values of throttle setting lies between points 
Pi and P 2 . The same notation is used to indicate the corresponding points on 
two of the trajectories in figure 6, 


Simplification of Minimum Fuel Profiles 

The minimiom fuel trajectories consist of acceleration at level flight to 
a speed which is nearly independent of the range, climb at maximum power set- 
ting, a segment during which the power setting is reduced from maximum to 
idle, and a descent at idle (see figs. 6 through 8), In addition, if the 

range to be covered is larger than h51 km, a cruise segment at = 226 m/sec 

and Hq = 10 km will be present. 

It has been found that irrespective of a given range, the minimum- fuel 
flight profiles can be replaced by approximating trajectories that consist of 

only two segments: (l) a climb at maximum climb power setting and 152.5 m/sec 

equivalent airspeed (EAS) and (2) a descent portion at idle and 128 m/sec 
(EAS). The change from the climb to the descent portion is made at an energy 
level which gives the desired range. The total fuel consumption associated 
with the approximating trajectories is about 1.5 percent higher than that of 
the minimum- fuel flight profiles. 


Minimum Time Profiles 

Minimum time trajectories, which correspond to a = 0, are shown in fig- 
ures 6, 7 9 and 8 along with trajectories for the other a values. For this 
case, figure 8 shows that power is either at idle or at maximum. The switch 
from maximum to idle occurs in the decreasing energy regimen when the airspeed 

has decreased to Vg , the airspeed at the maximum energy point. During a 

max 

part of the descent the airspeed follows the design dive flight placard of an 
EAS of 238 m/sec. 
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The structure of minimuiii time flight profiles can also he derived from 
the observation that the function to be minimized at each energy level is 
[1 - (V/VEjj^^)]/|v(T-D)| where Vjj^ax " velocity associated vith 

the maximum energy of the particular trajectory. During descent, if V < ^Ej^ax 

the minimizing value of the power setting is at idle; on the other hand, if 
V > Ve 9 the minimizing value is the maximimi power setting. 


Minimum DOC Profiles 

By letting distance be a parameter, a spectrum of optimum trajectories 
from minimum time to minimum fuel may be obtained as a is varied from zero 
to 1. The time-fuel tradeoff generated by this spectrum is plotted in figure 
9 for several distances. It is characteristic that these tradeoff curves 
approach a horizontal slope near a = 1 and a vertical slope near a = 0. 

Given the relative dollar costs of time and fuel, one can use figure 9 to cal- 
culate the time-fuel tradeoff that minimizes the DOC for a given range. As an 
example, assume $132/tonne fuel cost and $5/min time cost, typical costs for 
operating this particular aircraft. This yields a cost ratio between fuel and 
time of C<p/Cp = 5/132 = 0.038 tonne/min. The minimum DOC operating point is 


then obtained where the line with slope 



Figure 9.- Time-fuel tradeoff for fi 
ranges; CTOL aircraft. 


-0.03o tonne/min is tangent to the 
fuel-time tradeoff curve for the 
desired range. For a range of 
400 km, this is point A in figure 9 
and is seen to require only about 
3 percent more fuel than the minimum- 
fuel performance point. 

It is also possible to relate 
the parameter a to the cost ratio 
using equation (5), but then care 
must be exercised that the same units 
for mass and time are used in equa- 
tions (4) and (5) and that no scale 
factors are later introduced in the 
algorithms for minimizing equation (4) . 
Such scale factors were used in the 
calculation and therefore the inter- 
mediate values of a shown in fig- 
ure 9 as well as those in earlier 
figures should not be correlated 
with C»p/C|p using equation (5). 
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Fixed-Range Optimum Flight Paths: STOL Aircraft 

Except for minor variations, the same methods vere used for the STOL air- 
craft as for the CTOL. Because the STOL model used is the result of a prelim- 
inary design study, it is important to note that the actual numerical values 
presented might be subject to considerable revision in the future. Neverthe- 
less, the results serve to illustrate the effects vhich are peculiar to an 
aircraft equipped for the use of powered lift. 

The method of keeping the thrust setting vithin the allowable limits dif- 
fered somewhat from the case of the CTOL aircraft. The turbine inlet tempera- 
ture Tin ’W’as limited in the optimization to lOOO’^K < Tin ^ lH50^K and the 
corrected turbine temperature, Tq = Tin/Q, was limited as shown by the "max 
Tc" boundary in figure 4(a), Regardless of which constraint is governing, the 
lower limit will be referred to as "flight-idle" and the upper limit as 
"maximum thrust," 


The optimization for this aircraft over most of the descent portion 
(except when a = O) was complicated by the existence of two local minima in 
the integrand of equation (l9). These local minima varied, from being widely 
separated to nearly coincident, depending on the energy level. The function 
minimization routine located both local minima, and then chose the controls 
corresponding to the smallest of the two as the optimum controls. 


Figure 10 shows the contours of cruise fuel efficiency versus true air- 
speed for different values of altitude using an aircraft weight of 90,700 kg. 
The minimum fuel cruise condition occurs at an altitude of 11 km and a speed 
of 200 m/sec. As in the case of the CTOL aircraft, this point is below the 


ceiling; however, the difference is only 



TRUE AIRSPEED, m/sec 


Figure 10.- Cruise fuel efficiency; 
STOL aircraft. 


about 100 m in this case. The minimum 
time cruise condition occurs near 
an altitude of 6 km and an airspeed 
of 230 m/sec, A tradeoff between 
time and fuel can be obtained by 
selecting different operating 
points along any of the altitude 
contours. At all altitudes, the 
curves are quite flat near the mini- 
mum fuel point, and a speed change 
of about 7 percent is required to 
produce a 1-percent change in fuel. 
The fuel consumption near the mini- 
mum time curve increases more 
rapidly with speed, and an increase 
of 5 percent in time produces about 
an 8 percent reduction in fuel, 
except at the highest altitudes. 
Also, along the minimum-time curve 
the fuel consumption increases 
rapidly with decreasing altitude 
while corresponding changes in 
speed are small. 
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A series of minimuni-fuel and miniin-uiii-time trajectories was calculated for 
various values of Emax including the minimum- fuel and minimum-time 

cruise energies. The initial and final energies were 1,5 km in all cases, and 
a lower altitude limit of 150 m was Imposed. The results for intermediate 
values of a were not determined for this aircraft. 


Minimum Fuel Profiles 

Figure 11 gives the altitude-range profiles for three minim-um-fuel tra- 
jectories having no cruise section. The maximum energy for the longest trajec- 
tory is the minimum-fuel cruise energy. Note that all of the trajectories fol- 
low the same ascent path until maximum energy is reached. The flight-path 
angle appears to be discontinuous between ascent and descent at the scale used 
in the figure, but the use of sufficiently small increments in energy shows 
that this is not the case. For short ranges, most of the distance is covered 
during the descent portion, but as the range, and hence the maximum energy 
increases, a greater portion of the distance is covered during ascent. This is 
because the energy rate and the corresponding flight-path angle approach zero 
slowly near the minimum- fuel cruise energy. The final portion of each trajec- 
tory consists of a steep descent followed by a short section of level flight 
at the low altitude limit. 

The speed-altitude profiles in figure 12 are for the minimum- fuel trajec- 
tories just discussed and two minimum-time trajectories. In the minimum- fuel 
case, the ascent and higher altitude portion of the descent segments are quite 
similar to the results presented earlier for the CTOL aircraft. In this 
region, maximum thrust is used for climb, and flight-idle is used for descent, 
except near the maximum energy, where a gradual transition is made to or from 
the cruise thrust. In the lower portion of the descent, the airspeeds become 
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Figure 11.- Altitude-range profiles, 
STOL aircraft. 


Figure 12.- Speed-altitude profiles, 
STOL aircraft. 
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very high and thrust settings greater than flight-idle are used to help attain 
the high speed, after vhich the thrust is again reduced. 

As the aircraft approaches the maximum speed on the descent trajectory 
(fig. 12 ), the magnitude of the flight-path angle and the magnitude of the 
acceleration hecome very large. For example, consider the case of the trajec- 
tory for the longest range ( 76 ^ km). A 500-m reduction in energy hetveen 
points A and B occurs in 2k sec and is accompanied hy a speed increase of ^3 
m/sec and an altitude change of nearly 1,6 km. During this interval the 
flight-path angle reaches -37°. Although the diving segments hecome shallower 
at lower values of Emax? trajectories appear to approach the constant 
energy transitions encountered in other studies (ref. 2), 

The steeply diving portions are prohahly caused by the high values of 
flight-idle thrust, and hence fuel flow to the engines providing powered lift 
for this STOL aircraft. At low altitudes, where the fuel consumption is high, 
fuel usage is minimized by descending at high speed and then decelerating to 
the final energy as rapidly as possible. 

The high speeds at low altitudes and the large decelerations and flight- 
path angles encountered in the descent are unacceptable for passenger aircraft. 
One possible constraint, which provides an acceptable suboptimum trajectory 
for approximating the minimum- fuel case, is illustrated by the dashed lines in 
figure 12. For ascent, the indicated airspeed was restricted to a maximum of 
130 m/sec (250 knots) and for descent the limit was chosen to fit the envelope 
of the upper portions of the descent segments. This limit provides fairly 
smooth transitions between the constrained and unconstrained portions. 


Minimum Time Profiles 

The minimum-time profiles in figure 12 are similar to those shown for the 
CTOL aircraft. The longest range trajectory (l064 km) just reaches the mini- 
mum time-cruise energy. The horizontal portions result from constraining the 
altitude to be at least 150 m; the initial and final energy levels were 1.5 km. 
Maximum thrust is required almost during the entire trajectory, except during 
the final horizontal segment for V < when it must be at flight-idle. The 
maximum speed during descent is limited by the drag increase with Mach number 
rather than by a control constraint. The magnitude of the flight-path angle 
never exceeds 


Fuel and Time Tradeoff of Trajectories 

The fuel and time used for various ranges for the minimum-time, minimum- 
fuel, and constrained-trajectories are plotted in figure 13, with the ranges 
indicated on the curves. The curves for intermediate a values were not 
obtained for the STOL aircraft. The effects of constraining airspeed in the 
minimum-fuel case are shown by the dashed line in figure 13, The constraint 
increases the cost in both fuel and time at each range. The increased costs 
occur almost entirely during descent and are attributed to the elimination of 
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TIME, min 

Figure 13.- Fuel and time require- 
ments for fixed ranges , STOL 
aircraft. 

descent required about 10 percent more 
minimum- fuel cruise and descent coveri; 


the high-speed segment during descent. 

A less severe speed constraint for 
descent (e.g., the one used for ascent) 
would provide a cost nearer the optimum 
if a satisfactory transition could be 
made to the optimum. 

Since the ascent portions of the 
optimum and constrained trajectories 
are almost identical, it is desirable 
to compare costs only during the 
descent portion. The lower airspeeds 
in the constrained trajectory resulted 
in lower magnitudes of flight-path 
angle and energy rate, which in turn 
produced a larger range and longer 
flight time for descent from a given 
maximum energy. For the case of the 
maximum energy equal to the optlmimi 
cruise energy, the ranges were equal- 
ized to 891 km by adding an appropriate 
cruise section to the optimum trajec- 
tory. In this case, the constrained 
fuel and 19 percent more time than the 
g the same range. 


CONCLUSIONS 


By using the aircraft’s total energy (kinetic plus potential) as the inde- 
pendent variable in the climb and descent segments, an efficient algorithm has 
been developed for computing fixed-range optimum trajectories for short-haul 
aircraft. The good computational efficiency of the algorithm makes it attrac- 
tive for parametric studies as well as for possible on-board implementation. 

In this report, fuel cost and flight range were selected as parameters and used 
to compute optimum time-fuel tradeoff curves for the in-service CTOL aircraft. 
Such curves can help an airline in determining the best operating strategy in 
periods of fluctuating fuel prices. The optimum trajectories for both the CTOL 
and STOL aircraft are characterized by the absence of a cruise segment unless 
the range of the mission exceeds a minimum value. For the CTOL aircraft, the 
optimum controls, consisting of airspeed and engine power setting, were smooth 
functions of the energy, thereby yielding fly able, though nonstandard, trajec- 
tories. However, additional work is required to define the on-board displays, 
computers, and autopilot navigation system interfaces required for the pilot 
to be able to fly the optimum trajectories. For the STOL aircraft, the air- 
speed was nearly discontinuous with respect to energy at one energy level dur- 
ing the minimum fuel descent. Under the worst conditions, elimination of the 
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steep dive, required to fly such a descent, increased fuel consumption by 10 
percent and the time to fly by 19 percent compared to the optimum descent 
trajectory. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., 9^035, June 3, 1975 
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APPENDIX 


PROPERTIES OF THE OPTIMUM CONTROLS NEAR MAXIMUM ENERGY 


In this appendix, we study the properties of the two functions (l8) and 
( 19 ) to be minimized at the maximum energy Ejj^ax* This will lead to conditions 
under which equations (21a) and (21b) hold. In the discussion that follows, 
thrust T, rather than throttle setting II, and velocity V are chosen as 
controls. The region of allowable controls at E = Ej^g^^ is denoted by Q. 

The cruise controls {T^,, obtained from equation (17/ lie in the interior 
of Q except for values of a near zero, where thrust or airspeed or both 
are at their maximum allowable values. In this analysis, we only consider the 
case where {T^, V^} are in the interior of Q. Standard subscript notation is 
used to identify partial derivation with respect ’to the controls T and V, 

As illustrated in figure 2, it can be shown that the gradients with 
respect to the controls of the functions P - ^opt^ V(T-D) are parallel at 

the cruise point, denoted by Q. This is expressed as 

^ “^V ^opt 

Defining the quantities inside the minimization operators of equations ( 18 ) 
and ( 19 ) as g: 


4 P - 

E 


Near Q, the function g can be represented approximately by a second-order 
Taylor series expression in the controls: 


g(T,V,E ) 
’ max 


P^(dT - D^dV) + (I P^^dT2 + P^^dT dV + I P^dv2) 
(V^ + dV)(dT - dV - I P^dV2) 


(A2) 


where equation (Al) was used. The quantities dT and dV denote small incre- 
ments in the controls at the cruise condition. 


Local Properties of s(T,V,Ej^g^^) Near the Cruise Point 

Equation (A2) can now be used to study properties of g(T,V,Ejj^g^ 3 ^) as we 
approach the cruise point Q along any linear direction given by 


dT = 6 dV 


(A3) 
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or along the direction given hy 


dV = 0 (aU) 

Using L'H&pital's rule with equation (A2), we obtain 


y 


"Q 


lim 

dT=B dV 
dV^O 




'P^/V^ if B ?! 

-(P^B!!f2P„B*P„)/V^D^ if 



(A5) 


The direction in the control space defined hy 3 = I>y- must he excluded since 
it is the direction of the line through Q that separates the control space 
into the two regions ^ £ > 0 and E < 0. 

By considering g(T^V,E^g^y ) ~ gQ, it can he checked that Q is a local 
minimim of climh (i.e.^ if we approach Q with T > D) 

and that it is also a local minimum of in the descent (i.e., if 

we approach Q with T < D) along any direction given hy equations (A3) or (aU) 
with 3 ^ Dy, if and only if 

A = + 2(P^V^-P^)3 + [V^P^+P^(2D^+V^D^)1 > 0 for all 3 (a6) 

Equation (a 6) can he equivalently expressed as 


P >0 
TT 


(p V - VfVP +P(2D + VD )1<0 

^ TV c T^ TT o'- c W ^T^ V c 


(AT) 


Moreover, the minimum of g('^ 9 ^ 9 ^ax^ Pt/^c* 


As a consequence, if equation (AT) is satisfied, then in a small region 
near Q 


mln{g(T,V,E )} + min{g(T,V,E )} = 

T>D T<D ^c 


‘T 


V 


= 0 


(A8) 


thus verifying equation (21a). 

If it had heen assumed that P(T,V,Ejj^g^x) = where ^SYC 

the specific fuel consumption, as in some earlier work (refs. 2 and 3 ) , then 
equation (AT) is not satisfied. In fact, in this case the thrust values that 
minimize the two integrands for any E ^ E^^^ are a>nd T^_. ^ , respectively. 
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Global Properties of gCT^V^E ) 

max 

We will now give additional conditions under which Q will be the global 
minimum of during the climb and be the global minimum of 

-gCTjVjEinax) during the descent. 

If we assume that all the third and higher partials of P and D with 
respect to T and V are globally (in the control region of interest) negli- 
gible, then equation (A2) can be used to study global properties of g (this 
assumption was found to be reasonable for the CTOL aircraft considered in this 
report ) , 

By considering g - Pt/Vq it can be checked that in general g - P^/^c 
has a second local minimum corresponding to the velocity increment given by 

dV = (A9) 

The corresponding thrust increment is given by 

dT = -AB/f^^ (AlO) 


where A is given by equation (a6). 

For Pt/^c "to *t)e the global minimum of g In Q requires that this 
second minimum be outside A sufficient condition for this to hold is that 

(PtpyVc - “ ^TT^c^^c^W ^ ^ (All) 

where 


if 

Vw ’ » 

(this is the case for the CTOL aircraft) or 

X = V 

max 

W < ° 

The conditions shown in equations (AT) and (All) were verified for sev- 
eral cruise conditions for the CTOL aircraft using the analytical expressions 
for the propulsion and aerodynamic systems of the aircraft described in the 
body of this report. 
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